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ABSTRACT 


The  relationship  between  the  bending  moment  and  the 
through-plane  shear  force  in  the  vicinity  of  a mechanical  fas- 
tener at  failure  was  determined.  Experiments  were  conducted 
on  4-inch  wide  flat  plate  aluminum  and  graphite-epoxy  composite 
speciments  that  modeled  portions  of  a wing  skin  along  a spar 
and  along  a rib.  The  composite  specimens  were  either  8-ply 
or  16-ply  balanced  layups  and  were  simply  supported  at  two 
opposing  edges  and  free  along  the  other  two  edges.  The  fas- 
teners were  pulled  normal  to  the  plates,  and  the  maximum  force 
at  failure  was  measured  for  specimen  lengths  varying  from  two 
to  six  inches  between  supports. 

The  aluminum  plates  failed  by  formation  of  a plastic  hinge 
across  their  width  and  showed  little  sensitivity  to  through- 
plane  shear.  The  8-ply  spar  specimens  cracked  across  their 
width  and  also  were  relatively  insensitive  to  through-plane 
shear.  However,  failures  of  the  rib  specimens  were  confined 
to  a region  near  the  fastener,  where  the  fastener  pulled 
through  the  plate,  and  showed  much  greater  sensitivity  to 
through-plane  shear. 

Two  analyses  were  made;  one  for  small  elastic  deflections 
of  a thin  orthotropic  plate,  and  another  for  a beam  in  the 
elastic  range.  A mesh  generator  for  a finite  element  model  of 
the  plate  around  the  fastener  was  also  developed  for  the 
computer  program  ADINA. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  recent  years  considerable  effort  has  been  devoted  to 
the  development  of  new  materials  for  the  purpose  of  improving 
aircraft  performance.  One  family  of  these  materials  is  fiber 
reinforced  composites.  Many  studies  have  shown  that  use  of 
advanced  composite  materials  in  aircraft  structures  can  result 
in  significant  weight  savings  due  to  their  low  density,  high 
modulus  character.  Furthermore,  the  anisotropic  nature  of 
fiber  reinforced  composites  challenges  the  aircraft  designer 
to  exploit  their  directional  nature  in  order  to  realize  even 
greater  savings. 

Graphite-epoxy  composite  materials  have  recently  been  speci- 
fied for  major  structural  elements  in  two  U.S.  Navy  and  Marine 
aircraft  currently  in  advanced  development  stages.  Specifi- 
cally, the  wing  skin,  leading  edge  extension,  trailing  edge 
flap,  rudder,  horizontal  stabilizer,  and  vertical  fin  are  to 
be  constructed  of  graphite-epoxy  in  the  F-18,  [1].  The  under- 
lying wing  substructure  is  to  be  constructed  of  metal  in  a 
traditional  spar-rib  arrangement.  The  second  aircraft,  the 
AV-8B,  is  the  advanced  version  of  the  AV-8  Harrier.  This 
aircraft  will  have  the  entire  wing  box,  consisting  of  the  skin 
and  the  torque  box  substructure,  constructed  of  graphite- 
epoxy,  [2] 
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Mechanical  fasteners  are  specified  to  attach  the  wing  skins 
to  the  underlying  wing  substructure  for  both  aircraft  to  allow 
removal  of  the  skins, [3] . 

Extensive  tests  under  various  conditions  have  been  con- 
ducted to  assess  the  risks  associated  with  graphite-epoxy 
materials.  For  example,  strength  degradation  due  to  moisture 
absorption,  high  temperature  environments,  galvanic  corrosion, 
and  fatigue  have  been  determined  for  the  service  environment 
of  both  aircraft,  [3] . These  tests  have  shown  that  design 
strain  levels  for  both  the  Harrier  and  the  F-18  wing  designs 
provide  adequate  safety  margins  for  the  detrimental  effects 
of  elevated  temperature  and  moisture  absorption  to  be  encoun- 
tered in  service. 

In  addition  to  the  ability  of  the  composite  structure  to 
withstand  design  loads  under  a wide  range  of  service  conditions, 
some  attention  has  been  given  to  the  vulnerability  of  air- 
craft composite  structures  in  a combat  environment.  Tests 
have  been  conducted  to  study  the  response  of  advanced  composites 
to  ballistic  impact  damage  and  the  concomitant  reduction  in 
strength.  A discussion  of  these  tests  is  given  in  Reference  [4] . 

Much  data  are  available  and  a great  deal  is  known  about 
the  behavior  of  metals  at  failure.  Theories  of  failure  have 
been  tested  time  and  again  and  failure  models  have  proven  va- 
luable in  predicting  stress  levels  that  lead  to  failure.  With 
the  advent  of  composite  materials,  attempts  to  model  the  behav- 
ior of  composites  and  the  failure  mechanisms  have  led  to  some 


r 


1.1  U'*P<P  — 


understanding  of  failure  modes.  However,  significant  gaps 
exist  in  the  literature.  Of  particular  interest  here  is  the 
failure  of  a composite  structure  at  a metal  fastener  in  a 
particular  mode  called  fastener  pull-through. 


B.  FASTENER  PULL-THROUGH 

Consider  a composite  plate  or  skin  fastened  to  the  under- 
lying structure  with  mechanical  fasteners,  as  in  the  F-18  and 
Harrier  wing  designs.  Primary  considerations  in  selecting  skin 
thicknesses  include  in-plane  shear,  tensile  strength  and  criti- 
cal buckling  stresses.  Joints  are  primarily  designed  to  trans- 
fer in-plane  tensile  and  compressive  loads  in  the  skin  to  the 
underlying  spar-stiffener-rib  structure. 

Three  different  failure  modes  at  a jo^nt  are  shown  in  Figure 
1.  These  are  bearing,  shear-out,  and  tensile  failure.  These 
failure  modes  are  discussed  in  detail  in  Reference  [5] . A 
fourth  type  jf  joint  failure  is  illustrated  in  Figure  2.  In 
this  failure  mode  the  fastener  is  pulled  through  the  plate  by 
a large  through-plane  shear  force.  Although  very  little  atten- 
tion has  been  given  to  this  type  of  failure  in  metals,  this 
failure  mechanism  may  be  significant  in  fibrous  composites 
under  certain  loading  conditions. 

For  example,  through-plane  shear  stress  exists  in  the 
skin  of  pressurized  fuel  tanks.  Strain  levels  are  usually 
small  under  normal  design  conditions,  and  the  through-plane 
shear  stresses  at  the  joint  are  insignificant.  However,  if  the 
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aircraft  takes  a hit  in  combat,  circumstances  may  arise  which 
subject  joints  around  the  fuel  tanks  to  unusually  high  through- 
plane  shear  stresses.  The  extent  to  which  the  joint  is  able  to 
carry  these  stresses  may  have  significant  impact  on  the  vulner- 
ability of  the  aircraft. 

C.  HYDRAULIC  RAM 

Consider  the  problem  of  an  integral  wing  fuel  tank  which  is 
partially  full  of  fluid  when  impacted  by  a high  speed  projectile. 
The  kinetic  energy  of  the  projectile  is  transferred  to  the  fluid 
in  the  tank  as  it  passes  through,  [6] . The  result  of  this  en- 
ergy transfer  has  been  characterized  by  three  phases.  The 
first  phase  is  the  propagation  of  a shock  through  the  fluid. 

The  second  is  the  formation  of  a vapor  filled  cavity  behind  the 
projectile.  The  third  phase  is  an  oscillatory  phase  during 
which  the  cavity  expands  and  contracts  as  the  system  returns 
to  an  equilibrium  state.  These  phenomena,  illustrated  schema- 
tically in  Figure  3,  have  been  referred  to  as  the  hydraulic 
ram  effect. 

The  intense  fluid  pressure  created  by  the  hydraulic  ram 
has  a two-fold  effect  on  the  tank  wall.  First,  the  pressure, 
acting  perpendicular  to  the  wall,  creates  a moment  distribu- 
tion through  the  plate.  Secondly,  the  through-plane  shear 
forces  that  arise  must  be  reacted  by  the  fasteners.  This 
situation  is  shown  in  Figure  4. 
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Reference  [7]  describes  hydraulic  ram  tests  that  were 
conducted  on  fluid  containing  aluminum  tanks  and  one  manifesta- 
tion of  the  failure  of  the  tank  walls  was  an  "unzipping"  of 
fasteners.  Although  the  fasteners  that  were  used  in  these 
tests  were  not  typical  aircraft  fasteners,  questions  were 
raised  about  such  connection  failures.  If  aluminum  plates 
and  rivets  failed  in  this  through-plane  mode,  then  it  may  be 
possible  that  composite  plates  will  fail  in  the  same  way. 
Although  composites  are  generally  pound-for-pound  stronger  than 
aluminum  in-plane,  they  may  be  considerably  weaker  in  the 
through-plane  direction.  In  addition  to  the  bending  moment 
and  through-plane  shear  induced  by  the  hydraulic  ram  loading, 
aircraft  skins  will  be  in  tension,  compression,  and  in-plane 
shear  due  to  maneuvering  when  the  hit  occurs.  This  combina- 
tion of  stresses  can  lead  to  a significantly  lower  failure  load. 

D.  OBJECTIVES 

Several  aluminum  and  graphite-epoxy  plates  have  been 
tested  to  determine  the  sensitivity  of  failure  to  the  com- 
bination of  through-plane  shear  and  bending  moment  due  to 
hydraulic  ram  loading.  Specifically,  the  objectives  were  to 
determine  the  force  required  to  pull  a fastener  through  the 
plates  for  a range  of  values  of  the  bending  moment  at  the 
fastener.  A failure  curve  relating  pull  force  and  bending 
moment  similar  to  a yield  curve  for  limit  analysis  was  sought. 
Since  reinforcing  fibers  do  not  run  in  the  through-plane 
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direction  in  composite  plates,  it  was  expected  that  through- 
plane  shear  strength  would  be  low  compared  to  aluminum  plates. 

Macroscopic  failure  mechanisms  in  both  materials  were 
observed.  In  order  to  verify  experimental  observations,  two 
analyses  were  made.  One  analysis  was  for  small,  elastic  de- 
flections to  determine  the  magnitude  of  the  errors  in  moment 
distribution  associated  with  treating  the  specimens  as  beams 
instead  of  plates.  Another  model  applies  to  beams  in  the 
plastic  range.  Additionally,  a mesh  generator  for  use  in 
the  elastic-plastic  finite  element  code  ADINA  was  developed. 
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II.  EXPERIMENTAL  APPROACH  AND  SET-UP 


A.  THE  MODEL 

From  Figure  4 it  is  apparent  that  the  hydraulic  ram 
loading  on  a plate  creates  a complex  internal  stress  dis- 
tribution throughout  the  plate  and  particularly  at  the 
fastener.  To  further  complicate  matters,  hydraulic  ram  is  a 
dynamic  problem  involving  large  pressures  applied  over  a 
short  time  interval.  In  order  to  facilitate  the  experimental 
and  analytical  approach,  some  simplifying  assumptions  were 
made.  First,  the  dynamic  effects  of  the  hydraulic  ram  were 
not  addressed,  although  it  is  known  that  they  are  of  great 
importance  in  inducing  catastrophic  failures  at  fuel  tank 
boundaries.  Secondly,  it  was  assumed  that  there  would  be  no 
in-plane  stresses  present.  Again,  it  was  recognized  that  these 
will  play  an  important  role  in  joint  failures.  These  assump- 
tions are  justified  here  on  the  basis  that  this  is  a first 
study  in  an  ongoing  program.  Future  studies  will  consider 
these  effects. 

Figures  5A  and  5B  illustrate  the  modeling  of  the  skin 
with  the  test  plates.  The  underlying  structure  may  be  envi- 
sioned as  either  a rib,  Figure  5A,  or  a spar.  Figure  5B.  The 
difference  has  importance  only  with  respect  to  the  composite 
plates.  For  example,  the  zero  degree  ply  direction  was  as- 
sumed to  run  parallel  to  the  main  wing  spars,  [1] . 


Figure  5C  shows  the  loading  set-up  to  simulate  the  com- 
bination of  bending  moment  and  shearing  force  that  exist  at 
the  fastener  due  to  the  hydraulic  ram  loading.  The  edges  at 
y **  ± b/2  were  assumed  to  be  free  and  those  at  x = 0 and 
X = a were  assumed  to  be  simply  supported.  The  free  condition 
was  assumed  because,  near  the  ribs  and  spars,  moments  due  to 
the  hydraulic  ram  loading  would  be  smaller  in  the  y direction 
compared  to  the  moments  in  the  x direction. 

B.  TEST  EQUIPMENT 

In  order  to  load  a specimen  as  suggested  by  Figure  5C, 
a devide  was  constructed  in  the  Naval  Postgraduate  School 
Machine  Shop.  The  simple  support  condition  was  provided  by 
triangular  knife  edges  fashioned  from  steel  bar  stock.  A 
means  was  provided  by  which  the  supports  could  be  positioned 
from  one  to  five  inches  from  the  center  of  the  specimen  in 
one-inch  increments.  By  moving  the  supports  outward,  the 
moment  at  the  fastener  becomes  larger  for  a given  pull  force. 
Thus,  supports  close  together  provide  a relatively  larger 
pull  force  and  smaller  moment  whereas  supports  farther  apart 
give  a relatively  smaller  pull  force  and  larger  moment.  The 
supports  were  attached  to  a rigid  steel  base  to  facilitate 
attachment  of  the  device  to  the  moveable  head  of  a Riehle 
200,000  lb.  testing  machine.  The  device  is  shown  in  Figure  6 
with  an  aluminum  plate  in  place  for  a test  pull . The  fasteners 
were  passed  through  a hole  drilled  in  the  specimens  and  at- 
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tached  to  a Clevis  device  which  in  turn  was  attached  to  a 
Baldwin  SR-4  5000  lb.  capacity  load  cell. 

The  load  cell  was  connected  to  a locally  constructed 
Wheatstone  bridge.  A power  supply  and  digital  voltmeter 
were  connected  to  the  bridge  circuit  in  the  usual  way  and  cali- 
brated to  read  an  integer  multiple  of  the  load  applied  at 
the  fastener.  Bridge  output  was  also  monitored  on  a strip- 
chart  recorder  calibrated  to  read  pounds  of  pull.  The  strip- 
chart  recorder  was  most  useful  in  providing  a time  history  of 
the  applied  loading  and  assisted  in  recording  significant 
events  which  occurred  during  the  loading.  Some  aluminum 
specimens  were  fitted  with  strain  gauges,  which  were  wired  to 
a Wheatstone  bridge  switching  device.  Gauge  output  was  moni- 
tored by  a second  digital  voltmeter  calibrated  to  read  micro- 
inches per  inch.  An  additional  channel  of  the  strip-chart 
recorder  was  available  to  monitor  the  strain  gauge  output 
selected  at  the  manual  switching  station.  Figure  7 shows  the 
entire  experimental  facility. 


C.  TEST  SPECIMENS 

1.  Aluminum  Specimens 

The  aluminum  specimens  were  machined  from  sheet  stock 
to  a four-inch  nominal  width.  Lengths  of  the  specimens  were 
varied  from  seven  to  fourteen  inches,  depending  on  the  in- 
tended location  of  the  supports.  A hole  was  machined  in  the 
center  of  the  specimen  to  provide  a nominal  two-thousandths 
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inch  interference  fit  for  the  fastener  as  specified  in 
Reference  [8] . 

The  fasteners  that  were  used  for  the  tests  are  the 
Hi-Tigue  Hiloks  that  have  been  specified  as  skin  fasteners 
for  portions  of  the  F-18.  These  HLT-328-8-18  fasteners 
are  often  specified  for  applications  which  require  high  fatigue 
resistance.  Their  unique  Hi-Tigue  feature,  however,  was  not 
germane  to  the  tests.  Figure  8 shows  the  dimensions  of  the 
fastener. 


2 . composite  Specimens 

The  composite  plate  specimens  were  fabricated  by  the 
author  from  Hercules  AS-3501-6  twelve- inch  pre-preg  tape. 

The  tape  was  a nominal  ten  mil  thickness  with  low  resin  con- 
tent. Five  16-inch  x 16-inch  plates  were  constructed  using 
either  8-ply  [0/±45/90]g  or  16-ply  [ 0/±45/02/±45/90] s balanced 
layups  and  were  prepared  in  accordance  with  Reference  [9] . 

The  composite  preparation  equipment  is  shown  in  Figure  9. 

Four-inch  wide  test  specimens  were  cut  from  the  16- 
inch  square  plates  after  post-curing.  Some  8-ply  specimens 
were  prepared  such  that  the  0-degree  ply  was  perpendicular  to 
their  axis  (the  spar  specimens) . The  remaining  8-ply  and  the 
16-ply  specimens  had  the  0-degree  ply  parallel  to  the  axis 
(the  rib  specimens) . The  plates  were  drilled  and  fitted  with 
the  same  fasteners  as  the  aluminum  specimens,  although  inter- 
ference fitting  is  not  specified  for  this  type  of  material. 
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Room  temperature  material  properties  from  Reference  [5] 
are  tabulated  along  with  computed  flexural  rigidities  in 
Table  I. 


III.  ANALYTICAL  AND  COMPUTER  MODELS 


A.  ANALYTICAL  PLATE  MODEL 

The  test  specimens  are  idealized  as  thin  plates  with 
two  edges  simply  supported  and  two  edges  free,  as  shown  in 
Figure  10.  The  deflections  and  moment  distribution  can  be 
approximated  by  assuming  the  load  to  be  uniformly  distributed 
over  the  shaded  rectangle  of  area  uv  . The  governing  equation 
for  a balanced,  symmetric  laminated  plate  is  given  by  [6] . 


D11w'xxxx+2lDl!+2D««)w'xxyy  <2> 

+D2  2W'yyyy  =<5<x'y) 

where  w is  the  plate  deflection,  q is  the  applied  pressure 

load  and  subscripts  denote  partial  derivatives.  D , D , D 

11  12  22 

and  D are  stiffness  coefficients.  In  the  special  case  of 
6 6 

isotropy  we  have, 


D = D;  D 
2 2 12 


uD;  D = D(l-u)/2 
6 6 


The  solution  to  Equation  (1)  for  the  portion  of  the  plate 

PTS.t  shown  in  Figure  10  is  given  by  [9]  . 

w = ? { (a„+A_cosh  (miry/a  ) + (3) 

■•a . . m m 


+Bm (mTr-Va  )sinh  (nvir  y/a  ) } sin  (mir.x/a) 
where , 

a_  = 4Pa4  (-1)  (m-D/2  sin 

it  m D uv 
1 1 

and  a is  the  span  of  the  plate. 


(4) 


Considering  the  unloaded  portion  of  the  plate  beyond 


line  ts  , assume  a deflection  surface  of  the  form, 


w1  =Z{A1mcoshtmry/a)  + B1  Onity /a)  sinh(imry/a)  (5) 

+ C^sinhfln  Try/h  ) + :ry/a  ) cosh0n.7ry/a  ) } sinfa  Trx/k  ) 

The  six  constants  A , B , A1  . . ...Dl  . must  satisfy  the 

m m m m 

free  edge  boundary  conditions  at  y = b/2,  and  continuity 
along  line  ts  . Applying  continuity  conditions  at  y ■ b/2 
gives 


w * w1 ; w , = w 1 , • w . -h1  , ; w , =w 

• ' 9 YL*  9 w * ' W ' ' ' 
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It  is  shown  in  Reference  [9]  that 
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m = a (y  sinh  2y  -cosh  2y _) 
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B -B*  = a/2  (cosh  2y) 
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Two  more  equations  are  available  from  the  boundary  conditions 
at  y=b/2  in  order  to  solve  for  the  six  constants.  Along  the 
edge  y=b/2  the  moment  and  the  shearing  force  must 


vanish.  Hence 
M 
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Substitution  of  Equation  (4)  into  Equation  (9)  gives, 


{Bim> 


4 Pa1*  (-1)  Cm-D/2 


m5ir5D  uvK 


(1-p  )a  cosha  - (p.,-3)sinhci  (p-l)a  sinha  -2  cosha 
v 2 m m L m i m m u 


(P2-l)sinham 


(1-Pi )cosham 


-C1  fl-p  Isinha-D*  ( (1-p  )amcosha  +2  sinha_) 
mv  m mv  v mm  m' 
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where 

k = (1-p^  0--P2)am  + [2(p2*1)  + (p^l)  Cp  2 ‘ 3)  ] coshamsinham 
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The  quantity  that  is  of  particular  interest  here  is  the 
distribution  of  M along  the  line  x = a/2.  A plot  of  Mx/M 
versus  2y/b  is  shown  in  Figure  11  for  the  aluminum  specimen, 
and  each  of  the  three  types  of  composite  specimens  for  various 
values  of  a,  where 
M = 


B.  PLASTIC  LIMIT  ANALYSIS  FOR  A BEAM 

If  the  aluminum  specimen  is  considered  as  a beam,  then 
the  total  bending  moment  is  given  by  Reference  [12]  as 


M = b 


za(x)  dz 
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Where  b is  the  width  and  t is  the  thickness  of  the  beam. 

Making  the  usual  assumptions  that  plane  sections  remain 

plane,  the  strain  at  any  point  is  given  by 

e = kz  (13) 

Where  k is  the  curvature  of  the  middle  surface,  z * 0.  For 

moments  below  a certain  critical  value,  M , all  stresses 

e 

are  elastic  and 

a Ee  =Ekz  (14) 

Substituting  into  Equation  (12)  gives 

M — -j-j — b Eict 3 (15) 

If  the  moment  is  increased  until  at  the  outer  fibers, 
z = + 1/2,  is  equal  to  the  yield  stress,  cy  , the  maximum 
elastic  moment  is  given  by 

Me"-Tbt2ay  (16) 

Since  only  the  outermost  fibers  of  the  beam  are  yielded,  the 

beam  will  continue  to  carry  additional  load  through  the  cen- 
tral fibers  until  they  too  are  yielded.  The  situation  is 
shown  schematically  in  Figure  12.  The  limiting  moment  is 
Case  (D)  of  Figure  12  where  the  entire  cross  section  is 
yielded  and, 

M » J_c  bt2  (17) 

o 4 7 

The  moment  in  the  beam  is  largest  at  the  center  of  the 
span  and  its  value  is  Pa/4.  Substituting  Equation  (11)  into 
Equation  (17) , the  maximum  pull  that  the  beam  can  carry,  Py, 
is  given  by  p b t2 


(18) 


C.  MESH  GENERATOR  FOR  A FINITE  ELEMENT  MODEL 


The  theory  of  matrix  structural  analysis  is  discussed  at 
length  in  Reference  [13] , and  a discussion  of  various  three- 
dimensional  elements  can  be  found  in  Reference  [14] . The 
studies  presented  in  Reference  [14]  by  Clough  demonstrate  that 
hexahedral  elements  with  nodes  on  the  sides  adequately  repre- 
sent the  bending  of  a simply  supported  plate.  Clough  used 
a single  layer  of  20  node  elements  and  achieved  the  exact 
deflections.  He  showed  that  more  complicated  elements,  such 
as  curved  tetrahedra,  could  perform  as  well,  but  that  the 
formulation  time  was  significantly  greater.  Further,  it  was 
shown  that,  while  two-point  Gauss  quadrature  rules  gave  exact 
results  for  rectangular  prisms,  three-  or  four-point  rules 
were  required  for  skewed  elements  in  order  to  achieve  the 
desired  accuracy. 

In  this  study,  16  node  isoparametric  elements  were  chosen 
to  model  specimens.  The  objective  of  the  finite  element 
effort  was  primarily  to  develop  a mesh  generator  which  could 
later  be  used  to  analyze  the  laminated  anisotropic  plates. 

The  finite  element  code  used  was  ADINA,  developed  by 
Bathe  [15]  in  1975.  ADINA  incorporates  a library  of  four 
basic  elements  and  allows  selection  of  any  combination  of  15 
material  models.  The  library  material  model  that  seemed 
best  suited  was  the  isotropic  elastic-plastic  model  with  strain 
hardening.  Tensile  tests  were  conducted  on  aluminum  coupons 
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of  three  different  alloys  in  order  to  accurately  provide 
proper  input  data  for  ADINA  . The  results  of  the  tensile 
tests  were  used  throughout  the  investigation  as  baseline 
material  properties,  and  the  stress  strain  curves  for  three 
different  aluminum  alloys  are  shown  in  Figures  13A  through 
13C.  The  bilinear  approximation  to  the  stress  strain  be- 
havior of  the  materials  is  indicated  by  the  broken  lines. 

Symmetry  of  the  test  specimens  allowed  some  simplifi- 
cation of  the  model.  The  analysis  includes  one-quarter  of  the 
plate.  The  finite  element  model  of  the  specimen  is  shown  in 
Figure  14.  The  loads  can  be  applied  as  concentrated  loads  at 
nodes  along  lines  (A)  and  (B)  of  Figure  14.  At  x * 0,  re- 
straints can  be  applied  such  that  no  displacement  in  the  x 
direction  is  permitted.  At  y - 0,  no  displacement  in  the  y 
direction  is  permitted.  At  x = a,  the  simple  support  condition 
dictates  that  no  displacement  in  the  z direction  is  permitted. 
The  only  other  constraints  involve  the  nodes  at  z = 0 and  on 
the  first  ring  of  elements  closest  to  the  fastener.  The  assump- 
tion is  made  that  the  fastener  is  essentially  rigid  with  respect 
to  the  aluminum.  Consequently,  no  displacement  in  either  the 
x or  y directions  is  permitted  for  these  nodes. 

A description  of  the  mesh  generator,  and  its  use  is  in- 
cluded in  Appendix  A.  Present  limitations  internal  to  the 
ADINA  source  code  preclude  the  possibility  of  obtaining  a 
successful  nonlinear  analysis  due  to  array  dimension. 


IV.  RESULTS  AND  DISCUSSION 


A.  ALUMINUM  SPECIMENS 


All  of  the  aluminum  specimens  tested  failed  by  forming  a 
plastic  hinge  across  the  center  of  the  span,  as  shown  in 
Figure  15.  None  of  the  fasteners  pulled  through,  and  there  was 
no  visual  indication  that  the  through-plane  shearing  stresses 
had  any  effect  on  the  failure  mode. 

Table  II  contains  the  ultimate  moment  and  pull  force 
based  upon  the  plastic  limit  analysis  and  the  actual  maximum 
moments,  Mf,  and  pull  forces,  Pf,  obtained  during  the  experi- 
ments. Note  that  due  to  the  fact  that  a » 4 inches,  Pf  and 
Mf  have  the  same  numerical  value.  The  ultimate  moments  ob- 
tained experimentally  are  seen  to  be  somewhat  higher  than 
those  predicted  by  the  plastic  limit  analysis  of  section  III. 
Part  of  the  difference  may  be  due  to  the  fact  that  the  theory 
applies  to  elastic-perf ectly  plastic  materials,  whereas  the 
materials  tested  (particularly  the  AL7075T6)  exhibited  some 
degree  of  strain  hardening.  Further,  the  limiting  moments 
predicted  by  the  theory  are  average  moments,  whereas  the 
moment  is  not  distributed  uniformly  across  the  beam,  as  shown 
in  Figure  11.  If  the  ultimate  pull  force  is  taken  to  be  that 
force  required  to  raise  the  moment  M along  the  Line  x = a./2 
to  the  ultimate  moment,  then  the  predicted  values  of  pull, 

Pf,  more  closely  resemble  those  obtained  experimentally. 

In  any  case,  the  tests  indicate  that  these  aluminum  plates  fail 
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at  the  joint  in  a plastic  hinge  mode.  This  hinge  effect  has 
previously  been  observed  in  hydraulic  ram  experiments,  [7] . 

Figures  16  and  17  show  a plot  of  the  stress  versus  applied 
load  for  two  of  the  locations  which  were  provided  with  strain 
gauges  on  the  6061  aluminum  specimen.  The  stress  is  normal- 
ized by  the  yield  stress  for  the  material,  and  the  load  is 
normalized  by  the  theoretical  Pf  for  the  material  and  support 
conditions.  The  figures  also  contain  the  theoretical  stress 
for  a beam  with  a concentrated  load  of  the  same  magnitude, 
and  for  the  flat  plate.  For  these  tests,  the  simple  supports 
were  4 inches  apart.  It  can  be  seen  that  at  the  point  x * 2, 

Y - 1.040,  the  beam  approximation  does  not  predict  stress  level 
accurately.  The  plate  theory  is  considerably  better.  Further 
away  from  the  centerline  of  the  plate,  at  x * 2.735,  y =.735, 
the  plate  theory  overestimates  the  experimentally  obtained 
stress  levels.  This  suggests  the  possibility  that  formulation 
of  the  plate  problem  may  not  be  accurate  near  the  fastener. 

The  assumption  that  the  fastener  head  behaves  like  a uniform 
load  acting  over  a rectangular  area  may  not  be  refined  enough 
to  accurately  predict  the  stresses  in  the  vicinity  of  the 
fastener. 

B.  COMPOSITE  SPECIMENS 

The  results  of  all  tests  on  the  composite  specimens  are 
presented  in  Table  III.  Mu^fc  was  determined  for  each  of  the 
five  plates  by  applying  a line  load  across  the  center  of  a 


four-inch  wide  by  seven-inch  long  specimen  cut  from  each  plate. 
The  line  load  was  applied  by  a steel  cylinder  attached  to 
the  Riehle  testing  machine,  as  shown  in  Figure  18A.  The  fas- 
tener passes  through  a hole  in  the  cylinder  and  the  hole  in 
the  plate  and  attaches  to  the  dlevis.  Thus,  pulling  the  fas- 
tener causes  the  cylinder  to  bear  against  the  plate  along  the 
axis  of  the  cylinder.  This  eliminates  the  concentration  of 
through-plane  shear  at  the  fastener.  Pul  was  determined 
for  each  16-inch  square  plate  by  pulling  a fastener  through  a 
two-inch  by  two-inch  specimen  restrained  by  a two-inch  by  six- 
inch  by  five-eighths-inch  thick  steel  plate  with  a circular 
hole  eighty-thousandths  of  an  inch  larger  than  the  radius  of 
the  fastener  head.  The  ?ult  set-up  is  shown  in  Figure  18B. 

Two  complete  sets  of  experiments  were  run  on  the  8-ply  rib 
specimens  in  order  to  establish  that  the  results  were  not  at- 
tributable to  defects  in  the  quality  of  the  plates  from  which 
the  specimens  were  cut.  The  mean  difference  was  less  than  7 
percent. 

Figure  19  is  a plot  of  the  experimental  data  in  nondimen- 
sional  form.  The  abscissa  is  the  moment,  M,  normalized  by 
Muit»  The  ordinate  is  the  pull  force  at  failure,  Pj,  normal- 
ized by  Puit*  The  extent  or  mode  of  failure  is  indicated  in 
the  figure. 


Figures  20A-20L  consist  of  photographs  of  the  8-ply  spar, 
8-ply  rib,  and  16-ply  rib  failed  specimens.  The  sequence  of 
photographs  for  each  group  of  specimens  is  from  small  moment 


and  mostly  through-plane  shear,  to  small  through-plane  shear 
and  mostly  moment  loading  conditions.  The  extent  of  the  failure 
zone  is  indicative  of  the  sensitivity  of  the  specimen  to  through- 
plane  shear.  Note  that  the  failure  zone  for  specimens  loaded 
with  relatively  small  moment  is  primarily  confined  to  a re- 
gion around  the  fastener  hole.  As  the  moment  was  increased 
(by  lengthening  the  span)  this  failure  zone  extended  in  the 
y direction  until  it  covered  the  entire  width  and  resembled  the 
failure  of  the  line  loaded  specimen. 

Figure  21A  is  a plot  of  the  data  from  a different  point  of 
view.  The  abscissa  in  this  plot  is  the  theoretical  elastic 
local  moment  at  the  fastener  head,  Mult,  normalized  by  Mult. 

Figure  21B  is  a plot  of  the  theoretical  elastic  moment  at  the 
edge  of  the  plate  at  the  center  cf  the  span,  normalized 

by  The  theoretical  values  were  obtained  from  Figures 

11A  and  11B  at  y = 0.4  inches.  In  both  cases  the  ordinate  is 
the  pull  force  at  Failure,  Pj,  normalized  by 

J 
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V.  CONCLUSIONS 

A.  ALUMINUM  SPECIMENS 

It  was  not  possible  to  pull  the  HL-328-8  fasteners  through 
the  aluminum  test  plates.  The  failure  of  the  aluminum  joints 
was  manifested  entirely  by  the  formation  of  a plastic  hinge 
and  the  resultant  inability  of  the  specimen  to  carry  further 
load. 

There  is  essentially  no  comparison  between  the  way  the 
aluminum  specimens  and  the  composite  specimens  failed  when 
subjected  to  the  combined  effects  of  moment  and  through-plane 
shear.  Failure  of  the  aluminum  plates  was  insensitive  to  the 
through-plane  shear  force  induced  by  the  head  of  the  fastener. 

The  assumption  that  the  plates  behaved  like  beams  when  sub- 
jected to  a force  on  the  fastener  does  not  appear  to  be  valid, 
with  the  exception  of  the  special  case  where  moment  (or  support 
spacing)  is  very  large.  It  was  shown  from  plate  theory  that 
the  effect  of  the  concentrated  load  upon  the  moment  distribution 
is  more  significant  where  support  spacings  are  small.  The  ap- 
proach  taken  here,  therefore,  requires  some  refinement  before 
it  can  accurately  predict  stress  levels.  However,  use  of  a 
finite  element  model  should  allow  designers  to  more  accurately 
predict  stress  levels  that  lead  to  failure.  A wide  variety 
of  loading  conditions  and  fastener  geometries  may  be  studied 
in  this  manner. 


B.  COMPOSITE  SPECIMENS  x 

It  is  evident  from  Figure  20  that  the  weaklv-oriented 
composite  spar  specimens,  52  and  53,  were  largely  insensitive 
to  through-plane  shear.  The  location  of  the  failure  points  on 
Figure  19  is  also  indicative  of  a lack  of  sensitivity  since 
these  specimens  failed  at  a load  that  was  very  near  the  ulti- 
mate moment.  However,  the  rib  specimen  showed  a significant 
sensitivity  to  through-plane  shear.  The  thickest  specimens 
showed  the  greatest  degree  of  sensitivity,  and  this  result  may 
be  cause  for  concern. 

The  manner  in  which  pull-through  failure  occurs  is  complex. 

It  is  suspected  that  the  ultimate  moment  is  attained  over  an 
area  near  the  fastener,  causing  failure  of  the  specimen  in 
that  vicinity.  However,  the  distribution  of  moments  over  the 
centerline  is  such  that  a portion  of  the  plate  outside  this 
failure  zone  maintains  sufficient  strength  to  carry  further 
load.  The  through-plane  shear  forces  must  still  be  carried  by 
the  composite  around  the  fastener.  Thus,  as  the  load  is  in- 
creased further,  the  fastener  pulls  through  the  specimen.  In 
those  specimens  which  were  loaded  with  large  moment  and  small 
pull-through,  the  distribution  of  moments  is  flatter  along  the 
centerline,  since  the  supports  are  further  apart.  That  is,  the 
effect  of  the  concentrated  load  induced  by  the  fastener  head 
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is  less  pronounced.  Conseuqnetly , once  local  failure  around 
the  fastener  occurs,  only  a small  increase  in  pull  force  is 
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required  to  increase  the  moment  to  its  critical  value  along 
the  entire  width  of  the  specimen.  Thicker  and  wider  speci- 
mens, in  which  the  ultimate  moment  is  larger,  may  show  greater 
sensitivity  to  through-plane  forces,  and  therefore  be  less 
likely  to  fail  across  their  entire  width.  It  would  be  in 
thick,  stiff  materials  where  "unzipping,"  if  it  is  to  occur, 
would  be  most  likely.  Current  wing  design  concepts  for  por- 
tions of  the  F-18  wing  call  for  thicknesses  from  .3  to  about 
0.7  inches.  Clearly,  as  thickness  gets  larger,  the  through- 
plane  shear  strength  gets  larger  and  there  is  a point  where  the 
actual  through-plane  shear  loads  can  be  carried  by  the  com- 
posite. 


VI.  RECOMMENDATIONS 


A.  ANALYTICAL  METHODS 

A mesh  generator  has  been  developed  to  model  the  simply 
supported-  simply  supported- free  - free  plate  loaded  over  a 
small  region  under  the  fastener  head.  This  has  been  done 
using  an  elementary  loading  condition  and  an  isotropic  material. 
A great  deal  of  understanding  could  be  gained  if  the  program 
were  modified  to  include  orthotropic  elements.  The  current 
program  allows  generation  of  up  to  16  layers  of  16-node  iso- 
tropic isoparametric  elements.  The  node  generation  and  ele- 
ment connectivity  would  remain  unchanged. 

A further  sophistication  could  be  realized  by  'treating  the 
nodes  near  the  fastener  as  fixed  and  applying  pressure  loads 
as  predicted  by  the  computer  program  Satans  [lfl.  Satans  has 
the  ability  to  predict  pressure  distributions  on  a plate  that 
arise  in  a fluid  as  a consequence  of  hydraulic  ram  from  a 
specific  threat  projectile.  Other  modifications  could  include: 

1.  The  ability  to  change  the  geometry  of  the  fastener  head. 

2.  The  ability  to  take  into  account  the  flexibility  of  the 
fastener  head. 

3.  Inclusion  of  in-plane  tensile,  compression,  and  shear  loads 
characteristic  of  maneuvering  flight. 

B.  EXPERIMENTAL  PROCEDURES 

No  data  were  obtained  at  the  lower  moment  ranges  because 
of  the  inability  to  move  the  support  arms  of  the  machine  closer 


together  than  two  inches.  This  design  defect  should  be  cor- 
rected by  affixing  the  support  arms  to  a base  which  permits  a 
full  range  of  motion,  and  does  not  restrict  selection  of  dis- 
tances between  arms  to  one  inch  increments. 

Investigation  of  the  sensitivity  of  composites  to  a more 
complicated  stress  state  which  includes  in-plane  stresses  is 
required.  Further  study  is  also  required  on  factors  which 
affect  sensitivity,  such  as: 

1.  Fastener  head  geometry  variations,  such  as  counter-sunk 
heads . 

2.  Fastener-hole  buffer  materials  and  seals,  such  as  O-rings. 

3.  Degree  of  interference  fit  or  tolerance  between  fastener 
and  hole. 

4.  Effects  of  high  cycle  fatigue. 

5.  Existing  delaminations  or  other  damage. 

6.  Effects  of  softener  strips. 


APPENDIX  A 


A.  GENERAL 

The  mesh  generator  is  extremely  simple  to  use.  Where  possi- 
ble, variable  names  are  exactly  the  same  as  in  Reference  [15] . 

In  addition  to  the  source  deck,  only  four  data  cards  are  re- 
quired to  generate  a data  deck  that  can  be  used  as  input  for 
ADINA.  By  simply  removing  the  load  cards,  the  deck  can  also  be 
used  with  the  preprocessor  program  PSAP1.  The  generator  is 
specific  and  therefore  limited  in  that  it  currently  has  the 
ability  to  generate  nodes  along  5 radial  lines  in  one  quadrant 
of  a plate  that  has  a hole  in  the  center.  Other  nodes  are  less 
restricted.  The  group  of  elements  closest  to  the  fastener  hole 
are  fixed  in  size  by  fastener  parameters  that  are  input.  The 
size  of  the  rest  of  the  elements  is  determined  by  the  number  of 
divisions  desired  along  the  line  that  lies  at  the  center  of  the 
span  of  the  plate.  The  capability  currently  exists  to  generate 


up  to  16  layers  of  16  node  isotropic  isoparametric  elements. 

However,  the  source  program  ADINA  requires  minor  internal  modi- 
fication to  accommodate  the  large  volume  of  storage  required  by 
a non-linear  analysis.  A program  is  currently  underway  at  the 
Naval  Postgraduate  School  to  provide  ADINA  with  this  increased 

i 

capability. 

B.  REQUIRED  INPUT  FOR  MESH  GENERATOR 

Input  data  for  the  generator  requires  only  four  cards. 


I 


Input  Card 

1 Format  (1615) 

COLUMN 

VARIABLE 

1 - 

5 

IDOF(l) 

6 - 

10 

I DOF (2) 

11  - 

15 

IDOF (3) 

16  - 

20 

IDOF (4) 

21  - 

25 

IDOF (5) 

26  - 

30 

IDOF (6) 

31  - 

35 

NEGNL 

36  - 

40 

MODEX 

41  - 

45 

NSTE 

46  - 

50 

IPRI 

51  - 

55 

NPB 

56  - 

60 

LAYERS 

61  - 

65 

NEGL 

66  - 

70 

NLOAD 

71  - 

75 

NLCUR 

76  - 

80 

NPTM 

Input  Card  2 FORMAT  (20A4) 


DESCRIPTION 


Number  of  layers  of  elements 
desired. 


1-80 


PTYPE 


ALPHANUMERIC  TITLE 

DESIRED  ON 

ADINA  AND  PSAPI  PRINTOUTS 


INPUT  CARD  3 FORMAT  (8 FI 0.0) 


1 

— 

10 

PROP (1) 

Young  1 s Modulus 

11 

- 

20 

PROP (2) 

Poisson's  Ratio 

21 

- 

30 

PROP (3) 

Simple  tension 

31 

- 

40 

PROP (4) 

Strain  Hardening  Modulus 

A Length  of 

41 

50 

half  span  between  simple 
supports 

51 

- 

60 

R1 

Radius  of  Fastener 

61 

- 

70 

PMAX 

Maximum  load  desired 

71 

80 

PINC 

Desired  load  increment  for 
load  data  cards.  Different 
from  DT  which  is  load  step 
internal  to  ADINA,  and  deals 
with  the  method  of  solution. 

INPUT  CARD 

4 FORMAT  (8F10 . 

0) 

1 

- 

10 

R 

Radius  of  Fastener  Head 

11 

- 

20 

T 

Thickness  of  specimen 

21 

- 

30 

DT 

Load step  increment 

Allow  to  default  (0.0) 

31 

- 

40 

TSTART 

Allow  to  default  (0.0) 

41 

— 

50 

FAC 

Allow  to  default  (0.0) 

39 


| 

i 

i 

i 


J 




nonnn  oo  ooooooo-g  o-  ui  onnn 


************* *****$****mesh  GENERATOR  PROGRAM  LISTING******* 


DIMENSION  ANGL(9) , ID!3,1500),X( 1500),Y( 1500 ),Z( 1500 ), 
*KEL(  32,  16 ) » KG!  32,  16)  ,N0D(  16,3  2»  16 ) , I PN0DE  ( 2 ,3 ) , I DOF  ( 6), 
*NPAR(20) « PROP (A ) , T IMV ( 20 ) , RV ( 20 ) , RV 1 ( 20 ) , NOOE ( 1500)  * 

READ! 2, 5)  ( ID0F<  I ) ,1=1 , 6 ) , NEGNL , MCDEX, NST  E , IPRI,NPB, 

*L AYERS, NEGL.NLOAD, NLC UR, NPT  M 
RE AD ( 2, 5 ) Nul  V » IDC,  I CON, I PR  I C , M AX ES 
FORMAT ( 1615) 

READt  2,6)  PTYPE 
FORMAT! 20AA) 

READ!  2,  7)  (PROP!  I ),I  = 1,A)  ,A  ,R  1 , PMAX,  PI  NC 
READ!  2,7)  R2,T,  DT ,T START  , FAC 
FORMAT! 8F10.0) 


ASSIGN  VALUES  TO  CONSTANTS  AND  SET  DEFAULTS 


NPR=2*NDI V+3 
NPRMID=(NPR+1 )/ 2 
DELR= ( A— R2 ) /FLOAT ( 2*NDIV ) 
N = 0 

NODE! 1 ) =00000 
NN=5*NPR+A*NPRMID 
NEL=A*1  NDIV+1) 

NDIV1=NDI V+l 
P I =4. 0* ATAN ! 1 • 0) 

ANGL ( 1 )=0.0 
ANGL ( 2 ) =ATAN! .25 ) 

ANGL ( 3 ) =AT  AN! . 5 ) 

ANGL ( A ) =AT  AN! • 75 ) 

ANGL! 5) =PI/A.O 

ANGL ( 6) =AT  AN! A. 0/3.0) 

ANGL ( 7 ) = AT  AN! 2 • 0 ) 

ANGL! 8) =ATAN( A. 0) 

ANGL  (9  )=P  1/2.0 
LAYRS=L  AY  ERS+l 
IPNODE! 1,1)=1 
IPNODE! 2, 1 ) =NN 
I F { NPB.  EQ.2)  GO  TO  8 
IPNODE!  1,  3)=NN  + 1 
I PNOOE! 2,3) =NN*LAY  ERS 
IPNODE! 1,2) =NN*LA  YERS  + 1 
IPNODE!  2,2)  =NN*LAYRS 
I F ( NLOAD. EQ.O ) NL0AD=23 
IF ( NLCUR . EQ  . 0)  NLCUR=2 
IF(NPTM.EQ.O)  NPTM=NSTE 
MTYP=i 
I ELD=0 
I ELX=0 
I ST=0 
1 0 IRN=  1 
NPTS=NP  TM 
NTF1=  1 
NT  F2  = 2 


GENERATE  THE  NODAL  MESH,  LAYER  BY  LAYER 
EACH  ELEMENT  WILL  HAVE  16  NODES,  8 ON  THE 


40 


Y-Z  UPPER  AND  8 ON  THE  Y-Z  LOWER  SURFACE. 

NO  NODES  ARE  GENERATED  THROUGH  THE  THICKNESS 
OF  THE  ELEMENT. 


DO  50  K=1  f LAYRS 

XN=T*<1  .0-FL0AT( K- 1 ) /FLOAT ( LAYER S) ) 

DO  40  J=1 ,9 
THETA=ANGL( J) 

I F ( MOD ( J , 2 ) .E3 . 0 ) GO  TO  21 
DO  31  1=  1 * NPR 
N=N+ 1 

IF(N.NE.l)  NODE ( N )=NQDE ( N- 1 ) +1 
IF(N.EQ.l)  NOD  E ( 1 ) = 1 
ID( L f N) =0 
I D(  2»  N ) =0 
I D < 3»  N) =0 

IF(J.EQ.l)ID(3tNl«I 

IF(J.EQ.9)ID(2»N)=1 

I F( J.GE .5. AND. I.EQ.NPR) ID(1  ,N)  *1 

IFIK.NE. LAYRS)  GO  TO  15 

IF(I.EQ.l)  I D( 2 ? N ) = 1 

IF( I.EQ.1)ID(3»N)*1 

X ( N)=XN 

IFU.NE.lt  GO  TO  16 

Y (N  ) = R l *COS (THETA ) 

Z(N)=Rl*SIN(THETA) 

GO  TO  30 

I F ( I.NE.2)  GO  TO  17 
Y ( N)  = ( Rl+(R2-Rl)/2.0)*C0S(THETA) 

Z (N)*(RH-(R2-R1 ) /2.0)  *SIN(  THETA) 

GG  TO  30 

IFU.NE.3)  GO  TO  18 
Y(N)=R2*C0S(THETA) 

Z(N)  = R2*SIN(THETA) 

GO  TO  30 

I FC  I.EQ.NPR)  GO  TO  19 

Y(N)=(R2+FLCAT( 1-3 )*DELR)*COS( THETA) 

Z (N )*( R2+FL0AT ( I - 3)*DELR) *SIN( THETA) 

GO  TO  30 

I F ( J. GT . 5 ) GO  TO  20 
Y (N)*A 

Z (N) =A*  TAN ( TH6T  A ) 

GO  TO  30 

Y (N)=A*CQTAN( THETA ) 

Z(N)=A 

IF(J.EQ.9)Y(N)=0.0 

CONT  INUE 

GC  TO  40 

DO  33  1=1 »NPRMID 

N=N+1 

NDINC=2 

IF(I.EQ.l)NDINC*l 
NODEIN) =NODE (N-l 1+1 
I D ( If  N ) =0 
ID(  2f N)=0 
I D(  3,N) =0 

IF  (J.GE  .5. AND.  I .EQ.NPRMID)ID(  L * N ) = 1 
IF(K.EQ.LAYRS.AND.I.EQ.l  ) 10  < 2,  N )=  1 
I D ( j,N) *ID( 2fN) 

X ( N I = XN 

IFII.NE.lt GC  TO  22 
Y(N)=Rl*COS(  THETA) 

Z ( N ) *Rl *S  IN (THET  A ) 

GO  TO  32 

I F ( I ,NE • 2 ) GC  TO  23 

Y ( N)  = R2*CQS  (THETA) 

Z(N)=R2*SIN( THETA) 

GO  TO  32 

I F ( I.EQ.NPRMID)  GO  TO  24 
Y(N)*(R2+FL0AT( I -2) *2. 0*DELR ) *COS ( THETA) 
Z ( N ) * ( R2+FLCAT ( I - 2 ) *2.0*DELR ) *S IN( THETA  > 
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24 


55 


56 


56 


59 


57 

99 

C 

C_, 

c 

c 

c 

c 

c 


GO  TO  3 2 

IFIJ.GT.5IGO  TO  25 
Y I N) *A 

Z(N)=A*TANI  THETA) 

GO  TO  32 

Y(N)=A*COTAN(THETA) 
Z IN  )*A 

IF ( J.EQ.9)Y(N)=0.0 
CONTINUE 
CONTINUE 
CONTINUE 


GENERATE  ELEMENT  CONNECTIVITY.  THIS  IS 
DONE  BY  SETTING  UP  THE  FIRST  ELEMENT  IN 
EACH  LAYER,  AND  REFERENCING  THE  REST 
OF  THE  ELEMENTS  IN  THAT  LAYER  TO  IT.  THE 
TASK  IS  SIMPLIFIED  TO  SOME  EXTENT  BECAUSE 
ADINA  HAS  PROVISIONS  FOR  ELEMENT  GENERATION. 


DC  99  1*1,  LAYERS 
KEL(  1 » I ) *( I -1 )*1  6+1 
KG( 1, I )*0 

NOD(l»l*I)=(I— 1) *NN+1 

NODI  2, 1 ,1 )=NOO< 1 , 1 ,1 ) +2 

NOD (3,1 ,1 )=NOD( 2 » 1, I ) +NPR+NPRM ID 

NOD (4,1,1) =NOD ( 3 , I , I ) -2 

DO  55  L = 1 » 4 

NOOIL+4,1,1 J=NOD(L,l, I)+NN 
CONTINUE 

NCD(9,1,I )*NOO(  1 , 1 , 1 ) +1 
NOD( 10,1,1) =NOD( 2,1,1 J+NPR—l 
NODI  11, l,I)*N0D(9.i.I J+NPR+NPRMID 
N0D( 12,1,1) *N0D( 10,1, I l-l 
DO  56  L*1 ,4 

NODIL+12, l,I)*NOD(L+8,l,I J+NN 

CONTINUE 

KOUNT  * 1 

00  57  J*2,NEL 

KCUNT=KOUNT+l 

KG(  J , 1 ) *0 

KELI J,II*KEL(J-l,I )♦! 

NOOI 1 , J ,1 J *NOD( 2 , J-l , I) 

IFIKOUNT.EQ.INDI  V+2) ) NOD( 1 , J , I ) =NOD( 4, J-ND I VI , I ) 

N00( 2 , J , I ) *NOD(  1 , J , I ) >2 

NOD( 3, J , I ) *NOD(  2 , J , I ) +NPR+NPRMI D 

NODI  4, J, I )*N00(3,J,I )-2 

DO  58  L«l,4 

NOO I L+4 » J » I ) *NOO I L,J,I)+NN 
CONTINUE 

NOD 1 9 , J , I ) * NOD (1,J,I)+1 
NOOI 10, J, I ) *NOD( 10, J-l, I )+l 

IF(K0UNT.EQ.(NDIV+2l)N0D( 10, J, I )*NOD( 10,J-NDIV1,I )+NPR 


NOOI 11, J, I ) *NOD( 3 * J, 1 1-1 
NOOI 12, J, I )*NOD( 10,J, II-l 
DO  59  L*1  ,4 

NODI  12+L,  J,  I)*N00(8+L,  J,I  M-NN 
CONTINUE 

IF(K0UNT.EQ.(N0IV+2) )K0UNT=1 
MAXE  L*KEL I J , I ) 

CONTINUE 

CONTINUE 


+NPRMID 


OIRECT  OUTPUT  TO  THE  APPROPRIATE  DEVICE. FOP 
DATA  CHECK  FUNCTIONS  THIS  WOULD  BE  THE  OFFLINE 
PRINTER • WHEN  THE  DATA  IS  SATISFACTORY,  OUTPUT 
BE  DIRECTED  TO  THE  OFFLINE  PUNCH.  IT  IS  THEN  RE 
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.(u  min  !•  iini'il  iii  . 


c 

c 


c 

160 

170 

180 

190 


200 

210 


220 


230 

240 


TO  BE  USED  WITH  THE  ADINA  JOB  CONTROL  DECK.  TO 
USE  IT  WITH  PS AP 1 SIMPLY  REMOVE  THE  LOAD  CARDS. 


WRI  TE  ( 7,160)MAXEL ,PTYPE 

FORMAT ( 1X.I5.2X. ’ ELEMENT  PLATE’ , 15A4) 

WRITE! 7,1 70) N,1  I DOF! I) , 1*1, 61 , N EGL ,NEGNL , MOCEX, NSTE , 


♦ OT.TSTART  »IPRI 

FORMAT ( 15,611,14,315, 2F 10.0,  15) 


250 

260 

270 

C 

C 

C 

C 

C 

c 

c 


280 


WRITE! 7,180) 

FORMAT!///) 

WRITE(7,190)  NPB, IOC 
FORMAT! 1615) 

WRITE! 7, 190 )i 1 IPNODEi  I , J ) , I * 1, 2 ) , J*1 ,NPB ) 

00  210  1*1, N 

WRITE!  7.200)  NODE! I ) , i 10 ( J, I ) , J* 1 ,3) ,X 1 I ) , Y 1 1 ) , l I I ) 

FORMAT! IX, I4,1X, 14,215, 15X,  3E10.3) 

CONTINUE 

WRITE! 7, 190  INLOAD, NLCUR,NPTM 

WRITE(7, 190) ICON, IPRIC 

DO  220  1=1,20 

NPAR! I )=0 

CONTINUE 

NPAR! 1 ) =3 

NPAR i 2 ) *MAX  EL 

NPAR! 3 ) *i 

NPAR!  7)*16 

NPAR(10)=3 

NPAR! ll)=3 

NPAR!  151=8 

NPAR! 16 )= 1 

WRITE!  7,230)! NPAR! I)  ,1=1,20) 

FORMAT! 2014) 

WRITE! 7 ,190)MTYP 

WRITE i 7 , 240 ) 1 PROP! I ), 1*1 ,4) 

FORMAT ( 8E10 .3 ) 

DO  270  1*1, LAYERS 
IPS=0 

IFU.NE.l  . OR.  I. NE. LAYERS)  IPS=1 

DO  260  J=l , NEL  , . . 

WRITE! 7,250)  KELi J , I ) , I ELD , I ELX , I PS, MTYP, MAXES, I ST , KG i J , I ) 
WRITE!7,190)(N0D! K,J,I),K=1  ,8) 

WRITE!  7,1 90 ) < NOD  ( K,J  ,1 ) ,K=9 ,16) 

FORMAT  18  15, El  0.3 ) 

CONTINUE 

CONTINUE 


LOAO  CARO  GENERATION  AND  OUTPUT  FOLLOWS 


DO  280  1=1, NSTE 
T I MV  ( I ) = I — I 

RV!I)*PINC/80.0*FL0AT! I) 

IFiRV!  I).GT.!PMAX/80.0))RV!  I)=PMAX/80.0 

RV1!I)*RV ( I ) / 2 . 0 

CONTINUE 

WRITE! 7,190  )NTF l , NPTS 
WRITE!  7,2  40 )( TIM V! I) ,RV!I  ),  1 = 1,  NPTS) 
WRITE! 7,190)NTF2, NPT  S 
WRITE(7,240)(TIMV1 I ) , RV 1 ( I) , 1*1, NPTS) 

DO  300  1 = 1,9 
NCUR= l 

IF! I.EQ.l .OR. I .EQ.9)  NCUR=2 
IF1M0D!  1,2)  .EQ.O)  GO  TO  290 
NLDl=LAYERS*NN  + l+  l I- I )*(  NPR+NPRMID1/2 
NLD2=NLDl+l 
NLO  3=NLD2+ 1 

WRITE! 7,3 10 )NLDl, 10! RN,NCUR , FAC 
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290 

300 

310 


WRITEI7.3  ^ . 

WRITE(7.310 )NLD3» I DI RN , NCUR , FAC 
GC  TO  300 
NLDl-NLDl+NPR 
NLD2*NL01 >1 

WRITE! 7.3101NLDI , 1 01 RN, NCUR. FAC 
WRITE! 7.310INLD2, IDIRN.NCUR  ,FAC 
CONTINUE 

FORMAT!  3I5.F10.0) 

STOP 


I0)NL02  , IDI RN. NCUR . FAC 
10 JNLD3, IDIRN.NCUR , FAC 


E 

l l 

* 

21.0  X 

10s 

psi 

E 

2 2 

» 

1.7  x 

10s 

psi 

G 

1 2 

= 

. 65x 

10s 

psi 

V 

1 2 

= 

.21 

V 

2 1 

* 

.017 

8-ply  spar 

8-ply  rib 

16-ply  rib 

Dn 

169.3 

634.4 

4815 

D 1 2 

99.7 

89.7 

914.1 

D 

2 2 

634.4 

169.3 

1379 

D 

6 6 

112.1 

112.1 

1014 

D 

1 6 

0 

0 

0 

D 

2 6 

0 

0 

0 

High  Strength  Graphite-Epoxy  material  properties 


TABLE  I 


j 


9 


Plate#  Specimen#  a Pull  Load  Moment  at  Mult  Pult  Remar 

at  failure,  failure. 
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Figure  18 


Figure  19 


small  moment,  16-ply,  P =1255  lbs. 

Pf  = 2925  lbs. 


(K)  a=4",  16-ply  rib  specimen 
P,  =1140  lbs. 


( L)  a=  5",  16-phy  rib  specimen 
P,  =1085  lbs. 


(M)  a=6',  16-ply  rib  specimen,  (N)  a=4",  16-ply  rib  specimen 
P.  =^85  lbs.  line  looded,  P.  =1780  lbs. 
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